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ABSTRACT 


Computer  simulation  studies  of  two  frequency  domain  adaptive  beamforming  algo¬ 
rithms  are  presented.  The  two  algorithms  are  the  frequency  domain  adaptive 
beamforming  modified  least-mean-square  algorithm,  and  the  proposed  new  frequency 
domain  adaptive  beamforming  algorithm  based  on  modified  adaptive  linear  prediction- 
error  filtering.  The  simulation  studies  were  conducted  to  determine  the  multiple 
broadband  target  localization  capability  and  the  full  angular  coverage  capability  of  the 
two  algorithms.  The  number  of  iterations  that  the  adaptive  algorithms  took  to  reach  a 
minimum  estimation  error  was  determined.  The  algorithms  were  evaluated  at  several 
signal-to-noise  ratios.  Finally,  using  the  results  of  the  simulation  studies  a  comparison 
between  the  two  algorithms  is  performed. 
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I.  INTRODUCTION 


One  of  the  major  research  items  in  underwater  acoustic  signal  processing  is  to  de¬ 
velop  new  sonar  signal  processing  algorithms  capable  of  quickly  and  accurately  solving 
target  localization  problems.  Present  technology  dictates  that  several  lines  of  bearing  to 
a  target  must  be  obtained  before  a  sonar  fire  control  solution  can  be  computed.  Ob¬ 
taining  these  lines  of  bearing  is  a  time  consuming  and  often  dangerous  task  due  to  the 
increased  probability  of  counterdetection  and,  as  a  result,  evasive  maneuvering  and  de¬ 
fensive  action  on  the  part  of  the  target.  A  sonar  system  capable  of  providing  timely, 
accurate  target  localization  while  minimizing  own  ship  maneuvering  would  result  in 
longer  firing  ranges  and.  therefore,  a  reduction  of  the  threat  to  one's  own  ship. 

Frequency  domain  beamforming  is  accomplished  by  applying  appropriate  phase 
shifts  at  the  sensor  outputs  of  an  array  to  account  for  the  relative  propagation  delays 
of  a  signal  from  a  particular  direction.  The  phase-shifted  signals  from  all  sensors  are  then 
added  together  coherently  to  realize  the  full  array  gain.  Discrete  Fourier  Transform 
(DFT)  beamforming  is  the  usual  method  of  determining  the  direction  of  arrival  of  a 
plane  wave  signal.  A  discrete  number  of  direction  cosine  bins  are  formed  and  each  bin 
corresponds  to  a  discrete  direction.  If  the  number  of  direction  cosine  bins  is  large,  very 
fine  spatial  resolution  can  be  obtained.  The  phase  shifts  needed  to  cancel  the  relative 
propagation  delays  can  be  determined  adaptively. 

The  research  performed  in  this  thesis  is  a  continuation  of  the  work  of  Ziomek  and 
Behrle  [Ref.  1]  and  fully  evaluates  the  capabilities  of  the  complex,  modified  least-mean 
square,  frequency  domain,  adaptive  beamforming  algorithm  that  they  developed  via 
computer  simulation  studies.  Also,  a  new  frequency  domain  adaptive  beamforming  al¬ 
gorithm  based  on  modified  adaptive  linear  prediction-error  filtering  is  proposed,  which 
promises  better  resolution  in  low  signal-to-noise  ratio  environments  especially  when 
multiple  targets  share  common  spectral  lines  [Ref.  2.3].  The  computer  simulation  studies 
were  designed  to  test  the  two  algorithms  specifically  for  their  multiple  broadband  target 
localization  capability,  their  full  angular  coverage  capability,  and  their  angular  resol¬ 
ution  as  a  function  of  the  input  signal-to-noise  ratio  (SNR)  at  a  single  element  in  the 
array.  Using  lower  SNR  values  than  those  used  by  Ziomek  and  Behrle  [Ref.  1]  in  the 
evaluation  of  the  modified  L\1S  algorithm,  and  changing  some  assumptions,  such  as  the 
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number  of  array  elements,  sampling  rate,  and  number  of  iterations,  a  comparison  was 
performed  between  the  two  algorithms. 

Each  target  was  modeled  as  a  broadband  sound  source.  As  a  result,  the  frequency 
spectrum  of  the  output  signal  from  each  element  of  the  planar  array  contains  several 
frequency  components.  An  estimate  of  the  bearing  and  depression  angles  associated  with 
each  frequency  component  is  provided  as  a  result  of  processing  the  output  frequency 
spectrum  from  each  element  in  the  array  through  the  frequency  domain  adaptive 
beamforming  algorithms.  Therefore,  if  each  target  exhibits  at  least  one  unique  frequency 
component  (or  spectral  line),  then  all  targets  can  be  located. 

Full  angular  coverage  is  the  ability  to  localize  a  target  regardless  of  its  relative  po¬ 
sition  with  respect  to  the  array.  The  broadside  case  is  the  easiest  since  it  is  at  this  posi¬ 
tion  that  the  far-field  beam  pattern  beamwidth  is  its  narrowest.  The  endfire  case  is  the 
most  difficult  since  the  far-field  beam  pattern  beamwidth  is  the  broadest  at  this  point. 
The  full  angular  coverage  and  multiple  broadband  target  capabilities  were  tested  simul¬ 
taneously. 

Baseline  results  were  the  first  assembled.  They  are  defined  as  the  bearing  and  de¬ 
pression  angle  estimation  errors  (measured  in  degrees)  as  a  function  of  harmonic  num¬ 
ber,  and  the  number  of  iterations  of  the  two  algorithms,  for  the  "no  noise"  case. 
Identical  cases  were  then  run  using  additive,  zero  mean,  white,  gaussian  noise  to  corrupt 
the  output  signal  from  each  element  of  the  planar  array.  Average  bearing  and  depression 
angle  estimation  errors  were  plotted  as  a  function  of  the  input  SNR  at  a  single  element 
of  the  planar  array,  and  the  harmonic  number. 

Chapter  II  describes  the  theory  used  in  the  development  of  the  two  frequency  do¬ 
main  adaptive  beamforming  algorithms.  First  the  already  known  modified  complex 
least-mean-square  (LMS)  adaptive  algorithm  is  summarized.  Then  the  new  algorithm 
which  is  based  on  modified  complex  adaptive  linear  prediction-error  filtering  is  discussed. 
For  both  algorithms  the  construction  of  the  direction  cosine  estimates,  and  the  angle 
estimates  will  be  presented  in  detail. 

Chapter  III  contains  computer  simulation  results  and  an  explanation  of  these  re¬ 
sults.  The  results  are  presented  as  plots  of  the  average  estimation  errors  of  the  bearing 
and  depression  angles  at  three  distinct  signal-to-noise  ratio  levels  (-3dB.  -6dB,  -9dB)  for 
five  different  cases  for  each  of  the  two  algorithms.  The  different  test  cases  include  targets 
evaluated  at  broadside,  endfire,  random  placement,  and  targets  which  share  a  spectral 
line.  Conclusions  concerning  the  effect  of  harmonic  number,  sampling  rate,  number  of 
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array  elements,  number  of  iterations  used  in  the  two  algorithms,  and  SNR  are  made.  The 
Appendix  contains  tabular  numerical  data  for  all  results. 

Chapter  IV  will  present  final  conclusions  and  recommendations  for  future  research. 
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II.  ANALYSIS 


This  chapter  will  present  a  brief  analysis  of  the  two  frequency  domain  adaptive 
beamforming  algorithms  for  planar  arrays  which  are  evaluated  in  this  thesis.  The  goal 
of  both  algorithms  is  to  localize  multiple  broadband  targets  by  processing  the  output 
electrical  signals  from  a  planar  array  of  sensors  to  provide  estimates  of  direction  (both 
bearing  and  depression  angles)  and  frequency  content  of  the  multiple  acoustic  fields  in¬ 
cident  upon  the  array.  The  development  of  the  modified  complex  least-mean-square 
adaptive  beamforming  algorithm  for  planar  arrays  is  a  summary  of  the  analysis  section 
of  a  paper  by  Ziomek  and  Behrle  [Ref.  1],  However  the  modified  complex  adaptive  lin¬ 
ear  prediction-error  filter  beamforming  algorithm  for  planar  arrays  is  presented  for  the 
first  time.  As  a  theoretical  background  for  both  algorithms,  the  principles  of  adaptive 
beamforming  for  planar  arrays  in  the  frequency  domain  are  presented  in  the  next  sec¬ 
tion. 


A.  FREQUENCY  DOMAIN  ADAPTIVE  BEAMFORMING  FOR  PLANAR 
ARRAYS 

Consider  a  planar  array  of  M  x  X  (odd),  equally  spaced,  point  source  elements  lying 
in  the  XV  plane  where  M  and  X  are  the  total  odd  number  of  elements  in  the  X  and  Y 
directions,  respectively.  Let  the  random  output  electrical  signal  received  at  time  instant 
/  and  element  fm,n)  in  the  array  be  given  by 


r{l,m,n)  =>•(/,  m,n)  +  n{l,m,n),  l  =  -L' . 0 

m  =  (2.1) 

n  =  —  A' . 0,....Y'  , 


where  yfl,m,n)  and  n^l,m,rty  are  the  deterministic  signal  and  random  receiver  noise, 
respectively, 


L'  =  (L-  l)/2 


(2.2) 
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M'  =  (M-  1)12 


(2-3) 


A'  =  (-V  —  l)/2 


(2.4) 


and 


L  >  2K  +  1  (2.5) 

is  the  total  number  of  time  samples  that  must  be  taken  per  element  in  order  to  avoid 

aliasing  when  the  signal  y(l,m,n)  is  composed  of  K  harmonics  [Ref.  4].  If  we  designate 

the  length  of  the  data  record  (i.e.,  the  fundamental  period  )  recorded  at  each  element  in 

the  array  as  T0  seconds,  then  the  fundamental  frequency  (i.e.,  the  FFT  bin  spacing  )  is 

f0  —  -Jr-  Hz,  the  highest  frequency  component  that  can  be  calculated  in  y^l.m.n)  is  K/0 
*  0 

Hz.  and  the  output  received  signals  must  be  sampled  at  a  rate  of 


fs  =  LjT0  samples  sec,  (2.6) 

where  L  is  given  by  Eq.  (2.5)  [Ref.  4].  Taking  the  DFT  of  Eq.  (2.1)  with  respect  to  the 
time  index  /  yields  the  complex  frequency  domain  samples. 


R{q.m,n)  =  Y{q.m,ri)  4-  A  '(q,m,ri),  q  =  -L',...,0,...,Z/ 

m  =  (2.7) 

n  =  -.V' . 0.....V  , 

where  the  index  q  represents  the  harmonic  number.  If  the  acoustic  field  incident  upon 

l  •  ng  .  A 

the  planar  arrray  is  a  single  general  plane  wave  g(/  +  — ~c — )  propagating  in  the  ±  n0  di- 
reci  on  (see  Fig.  1),  where  nQ  is  a  unit  vector  and  g(t)  is  an  arbitrary  baseband  function, 
then  it  can  be  shown  that  the  output  electrical  signal  at  time  instant  /  and  element  (m.n) 
in  the  array  is  given  by  [Ref.  4  p.  160]: 


y{l,m,n )  =  g(lTs  + 


u0mdx  +  v0  ndv 
c 


(2.8) 
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,  \ 

A 

Figure  1.  General  plane-nave  field  g  (f  +  — p2-)  propagating  in  the  ±  n,  direction 
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The  corresponding  frequency  spectrum  is  given  by  [Ref.  2) 


,  ,  j2nqf0u0mdx  N  ,  _  j2nqf0v0ndy 

Y{q,nt,n)  =  Lcq  exp(  + - ~c - )  exp(  + - ~c - ) 


(2.9) 


where 


cq  =  Y{q,0fi)IL,  q  =  —  L' . 0 (2.10) 

are  the  complex  Fourier  series  coefficients  that  can  be  used  to  represent  g(t)  by  a  finite 

Fourier  series  with  K  harmonics  during  the  time  interval  |r|  ^  =  ~T~  . 

^  *  0 


u0  =  sin  60  cos  \l/0 


(2.11) 


and 


v0  =  sin  60  sin  \p0  (2. 1 2) 

are  the  dimensionless  direction  cosines  with  respect  to  the  X  and  V  axes,  respectively; 
90  and  t j/0  are  the  depression  and  bearing  angles,  respectively;  d,  and  dy  are  the  interele¬ 
ment  spacings  in  meters  in  the  X  and  Y  directions,  respectively;  and  c  is  the  speed  of 
sound  in  meters  per  second. 

Equation  (2.9)  is  the  basis  for  the  frequency  domain  adaptive  beamforming  algo¬ 
rithms  to  be  presented  in  this  thesis.  The  algorithms  presented  use  a  planar  array  to 
estimate  both  0O  and  \p0  for  each  harmonic  q  present  in  the  multiple  incindent  plane-wave 
fields  via  frequency  domain  adaptive  beamforming.  Therefore,  we  shall  use  the  following 
generalized  version  of  Eq.  (2.S) 


k 


uokmdx  +  vokndy 
c 


) 


(2.13) 


which  represents  the  output  electrical  signal  at  time  instant  /  and  element  (m,n)  in  the 
planar  array  due  to  several  incident  plane-wave  fields,  where  g*(r),  Uo*.  and  v0,  are  the  ar¬ 
bitrary  time  function  and  direction  cosines  associated  with  the  kth  sound  source  (target), 
respectively. 
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B.  FREQUENCY  DOMAIN  MODIFIED  COMPLEX  LMS  ADAPTIVE 
BEAMFORMING  ALGORITHM 


The  frequency  domain  adaptive  beamforming  algorithm  to  be  discussed  is  based 
on  processing  the  output  complex  frequency  domain  data  R(q,m,n)  from  all  M  x  X  ele¬ 
ments  in  the  planar  array.  We  begin  by  defining  the  complex  estimation  error  as 


«(?)  =  *(<?)-£(?) 


(2.14) 


where 


M'  S’ 

s(q)  =  "2jy"v'  ^  ^  \R(q,m,n)\  exp[  +jLR(q, 0,0)2  (2.15) 

m=-  \r  n=—Y 


is  the  reference  signal; 


\r  v 

Hq)  =  ' L  \Tx  Y  Y  <={q.m)d{q,n)R{q,m,n) 

_  •  V  (2.1 6) 

cT(q)R{q)d(q) 

LM.X 

is  the  estimate  of  s(q);  c(q.m)  and  d{q.n)  are  the  unit  magnitude  complex  weights  in  the 
X  and  Y  directions,  respectively;  o[q)  and  d{q)  are  the  M  x  1  and  N  x  1  complex  weight 
vectors  in  the  X  and  Y  directions,  respectively;  given  by 


Q(q)  =  Mq,  ~  m,-Aq-0) . c(q,M')f 


(2.17) 


and 


d{q)  =  [%,  -  X'),:,d(q,0) d(q, \')f 


(2.18) 
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respectively;  and 


m- 


R(q ,  -  A/'.  -  .V)~./?(<7,  -  M',0).../^,  -  A/',.V) 
R(qA-N')  ...  R(q, 0,0)  ...  /?(<?,(), A"') 

R(q,M\-V)  -  /?(<?, A/', 0)  ...  R(q,M',iV) 


(2.19) 


is  the  M  x  N  complex  data  matrix.  Next,  define  the  (M  +  N)  x  1  complex  weight  vector 
w(q)  as  follows: 


(2.20) 


Therefore. 


£(<?)  =  A  w{q) 


(2.21) 


where 


A=l  1  1  Q 

MxM  MxX 


(2.22) 


is  a  M  x  (\1  +  N)  matrix,  /  is  a  M  x  M  identity  matrix,  and  Q  is  a  M  x  N  null  matrix 
as  indicated;  and 


d{q)  =  R  w{q) 


(2.23) 


where 


R  =  l  Q  III 

A 'x\I  A’.rA’ 


(2.24) 


is  a  N  x  (M  +  N)  matrix,  Q  is  a  N  x  M  null  matrix,  and  /  is  a  N  x  N  identity  matrix  as 
indicated.  Substituting  Eqs.  (2.21)  and  (2.23)  into  Eq.  (2.16)  yields 
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(2.25) 


A/  ,  XT{q)Z(q)m 
S{q)  = - LMX - 

where 

Z(q)=ATRm  (2-26) 

is  a  (M  +  N)  x  (M  +  N)  complex  matrix.  The  complex  weight  vector  that  minimizes  the 
mean-square  error  £{|e(^)l2}  is  given  by  [Ref.  1] 

wi+1(?)  =  vv^)  +  2  W(?)[Z(?)  +  Zr(?)]V-  /  =0,1,2 .  (2.27) 

where 

ej{q)  =  s(q)  -  s,(q)  (2.28) 

is  the  estimation  error  after  the  ith  iteration. 


^(q)Z(qWq) 

LMS 


(2.29) 


is  the  estimate  of  s(q)  after  the  ith  iteration,  and,  the  step-size  parameter  q,  is  given  by 

+  '  =  0,1,2 .  (2.30) 

where  q0  is  constant  and  a}  and  <r2  are  the  signal  and  noise  power,  respectively,  at  the 
center  element  in  the  array.  After  each  iteration,  each  component  of  the  complex  weight 
vector  vv, _,(<?)  is  normalized  by  its  respective  magnitude  in  order  to  maintain  unit  mag¬ 
nitude. 

Once  the  complex  weight  vector  m,.,(q)  converges  to  a  steady-state  value  w^s(<7),  the 
steady-state  complex  weight  vectors  £ss(q)  and  dss(q)  can  be  obtained  from  Eqs.  (2.21) 
and  (2.23),  respectively.  The  estimates  of  60{q)  and  ij/0(q)  at  each  harmonic  q  are  given 
by  [Ref.  31 
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(2.31) 


d0(q)  =  sin-'[{[^s(?)]2  +  [tfS($)]2}1/2],  q  *  0 


and 


\j/0{q)  =  tan' 


aLS,  v 

vq  (q) 

*LS,  s 

uq  (q) 


qi*  0 


(2.32) 


where  u^{q)  and  v#5^)  represent  estimates  of  the  direction  cosines  obtained  by  using  a 
least-squares  fit  to  the  "unwrapped"  steady-state  phase  weights  Olss{q,m)  and  4>ss(<hn)i 
respectively.  Note  that  in  the  absence  of  noise 


6lss{q.m)  =  ± 


^qf0u0(q)mdx 

c 


q  =  -L',...rQ,...,L’ 
m  =  -M' . 0 . W 


(2.33) 


and 


4>Lss(q<n)  =  ± 


2^qfo^{q)ndx 

c 


q  =  -L',...,U,...,E' 
n  =  -V . 0 . X' 


(2.34) 


The  steady-state  phase  weights  need  to  be  "unwrapped",  that  is,  allowed  to  take  on  val¬ 
ues  outside  the  closed  interval  [  —  rr,  rc],  in  order  to  ensure  full  angular  coverage  (i.e., 
0  <  0Q(q)  <  ~  and  0  <  *p0(q)  <  2 n)  and  correct  depression  and  bearing  angle  estimates, 

A  A 

6a(q)  and  \Jj 0(q ),  respectively.  Finally,  substituting  Eqs.  (2.15)  and  (2.16)  into  Eq.  (2.14) 
yields  the  following  expression  for  the  steady-state  estimation  error: 


M'  .v 

ess(q )  =  l\[\  z  z  \R{q,m,ri)\  exp[  +jLR(q, 0,0)] 

m=— M'  n=—S" 

~  exP(  +jidlss(q^)  +  0ss(<7.«)  +  LR{q,m,n)']} 


(2.35) 


where 
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R{q.m,n)  =  \R{q,m,n)\  exp[  +jLR(q,m,n )] 


(2.36) 


Cssti,”')  =  exp[  =  exP[  +^55(<7.^)3  (2-37) 


and 


4ss(<7.«)  =  bss{q,n)  exp[  +j<bsS^<n)l  =  exp[  +j<bss{q,n) 3  (2-38) 


where  ass(q,m)=  1  and  bss(q,n)  =  1  are  real,  unit  magnitude,  amplitude  weights  and 
Bfs{q,m)  and  0j£( q,n)  are  real,  "wrapped"  phase  weights.  The  procedure  for  the  unwrap¬ 
ping  of  the  steady-state  phase  weights  in  order  to  obtain  least-squares  estimates  of  the 
direction  cosines  at  each  harmonic  is  described  by  Behrle  [Ref.  5]. 

C.  FREQUENCY  DOMAIN  MODIFIED  COMPLEX  ADAPTIVE  LINEAR 
PREDICTION-ERRROR  FILTER  BEAMFORMING  ALGORITHM 


Consider  a  planar  array  of  M  x  N  (odd),  equally  spaced,  point  source  elements  lying 
in  the  XY  plane  where  M  and  N  are  the  total  odd  number  of  elements  in  the  X  and  Y 
directions,  respectively.  Figure  2  shows  a  linear  prediction-error  filtering  structure  in  the 
frequency  domain  (i.e.,  after  the  operation  of  a  time  domain  DFT).  The  frequency  do¬ 
main  adaptive  algorithm  to  be  presented  is  based  on  processing  the  output  complex 
frequency  domain  data  R(q.m,n)  given  by  Eq.  (2-7).  Since  the  target  localization  prob¬ 
lem  requires  only  phase  informations  in  our  algorithm,  we  processes  phase  only  data. 

Let  the  one-step  forward  linear  prediction  errors  for  harmonic  number  q  in  the  X 
and  Y  directions,  be  given  by 


FpJ[t q,i,n)  =  sx(q,i,n)  -  sx{q,i,n),  n  =  -  A" . 0,...,.V 


(2.39) 


and 


Fp(qqn.i)  =  sy{q,m,f)  -  sy{q,m,i),  m  =  -  A/',...,0,...,M' 


(2.40) 
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1....  i 


i+P. ...  M  X 


h(q.i.n)  h£q,2,n)  h{q,3,n)  h£q,P.n) 


8(q.i.n) 


'  i  —  \  -m*  —  /  —  / 

i  -1 

J - *  r(q.i,n) 


s(q.i.n) 


Figure  2.  Frequency  domain  linear  prediction-error  filter  structure 


respectiv.  ly,  with 


iV_JLdL 

y 


„>  =  2Ld- 

M  2  . 


(2.41) 


(2-42) 


where  i  is  the  number  of  the  reference  element  in  the  direction  of  processing;  s,(  • )  and 
s,(  • )  are  the  reference  signal  and  the  estimate  of  the  reference  signal  in  the  X  direction, 
respectively;  and  sy(  • )  and  st(  • )  are  the  reference  signal  and  the  estimate  of  the  reference 
signal  in  the  Y  direction,  respectively.  Each  one  of  the  above  signals  is  defined  as  follows: 


sx(q,i,n)  =  R(q,i,n), 


sxiq,i,n)  =  y  Yhx{q,p,ri)R{q,i  +  p,n)  =  y  hl(q,n)dq,i,n), 
/=i 


(2.43) 


(2.44) 


13 


sy(q,m,i)  =  R(q.m,f), 


(2.45) 


and 

p 

A  1  V*  1  T 

sy(q,m,i)  =  —  2^hy(q,m,p)R{q,m,i  +  p)  =  yhy  ( q,m)z.{q,m,i )  (2.46) 

p=l 

where  h.(q,n)  and  h,(q,rrt)  are  the  P  x  1  complex  weight  vectors  in  the  X  and  Y  directions, 
respectively  ,  and  i{q,i.n)  and  z{q,m,i)  are  the  P  x  1  complex  frequency  domain  data  vec¬ 
tors  in  the  X  and  Y  directions,  respectively,  given  by 


hx(q.n)  =  [hJqA.n).  hx{q.2,n), . 

..,hx(q,P,n)']T, 

(2.47) 

r(q.i.n)  =  [R{q.i  +  1  ,n).R(q.i  4-  2.«),... 

. R(q.  ‘  +  Pm)Y, 

(2.4S) 

hv(q.ni)  =  [_hy{q.m.  1),  hy(q.m,2) . 

... ,hy{q.m.P)Y , 

(2.49) 

r(q.m.i)  =  [R(q,m,i  +  1  ),R(q,m,i  +  2),. 

. R{q.  m,i  +  P)~)T 

(2.50) 

where  P  is  the  order  of  the  linear  prediction-error  filter. 

Substituting  Eqs.  (2.43)  and  (2.44)  into  Eq.  (2-39)  and,  substituting  Eqs.  (2.45)  and 
(2.46)  into  Eq.  (2.40)  yields 

p 

FPt{q,i,n)  =  R{q-i,n)  -  -y  ^/i*(?,/’, «)£((?,/  +  p.n),  n  =  -  A',...,0 . .V',  (2.51) 

p=  i 

and 
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p 

Fp(q.m,i)  =  R{q.m,i )  -  ~  ^hY(q,m,p)R{q,m,i  +  p),  m  =  ~  A/',...,0 . A/'.  (2.52) 

p=i 

Since  we  are  working  with  phase  only  data,  Eqs.  (2.51)  and  (2.52)  become: 


and 


=  exPC  +jLR{q,i,n)~\ 

P 

~  "p  exp£  +JLhx(<bP’n)']  exp[  +jLR{q,i  +  />,«)],  (2.53) 

p=i 


Fpiq.m.i)  =  e\p[  +/Z.P(<7.m.';] 

p 

~  ~p  V|/il,(^,/n.ip)|  exp[  +/^/fv.(<7.wt,/>)]  exp[  +jLR(q.m.i  +  />)].  (2.54) 

/>=! 

m  =  -.\l' . 0,...,A/\ 

From  Eqs.  (2.53)  and  (2.54)  it  can  be  seen  that  if  each  component  of  the  complex 
weight  vectors  have  unit  magnitude,  then  the  forward  linear  prediction  error  become 
zero  when 


p 

exp[  +jLR(q.i.n)']  =  —  V  exp[  +jLhx{q,p,n)~\  exp[  +jLR{q,i  +  p,«)], 

p=\ 


n 


-.V  .  0  V 


(2.55) 


and 


exp[  +jLR{q.m,i)'] 


p 

1  ^ 


—  L  exp£  +JLh>^-mF)']  exp[  +jLR{q,m,i  +  />)], 

p=> 

rn  =  -A/' . 0 . A/'. 


(2.56) 
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Equations  (2.55)  and  (2.56)  will  be  valid  only  if  the  phase  weights  cophase  the  data  at 
the  output  of  each  element  in  the  prediction-error  filters. 

The  algorithm  used  for  the  optimum  estimation  of  the  complex  weight  vectors  is  a 
modified,  complex  least  mean-square  (LMS)  algorithm.  The  estimates  of  the  complex 
weight  vectors  in  the  X  and  Y  directions  after  the  k-th  iteration  in  X  and  V  directions 
is  given  by 


[4c(<M)L+i  *  Utx(<hn)lk  +  />)]*&:*(?, /.«)], 

n  =  —S' . 0.....Y'  (2.57) 

k  =  0,1, 2, 3 . 

and 


[Ay(^*w)]fc+i  =  [£v(<7-"0]k  +  2nklFPy{q,m,i)']k{r\q,m,i)'], 

m  =  -.)/' . 0 . M'  (2-58) 

k  =  0.1, 2.3 . 


respectively,  where  the  index  k  denotes  the  k-th  step  of  iteration,  Ft%{  • )  and.  F  ( .  )*  are 
the  forward  linear  prediction  errors  in  the  X  and  Y  directions,  respectively,  after  the  k-th 
iteration:  £*(  • )  is  the  complex  conjugate  of  the  data  vector;  and  the  step-size  parameter 
a*  is  given  by 


=  Mo  =  (4  +  alr\  k  -  0,1,2,...  (2.59) 

where  is  a  constant,  and  it  is  equal  to  the  inverse  of  the  sum  of  the  signal  and  noise 
power,  a *  and  a;  ,  respectively,  at  the  center  element  of  the  array.  After  each  iteration 
of  the  algorithm,  each  component  of  the  complex  weight  vectors  is  normalized  by  its 
respective  magnitude  in  order  to  maintain  unit  magnitude. 

Since  we  are  working  with  multiple  linear  arrays  to  obtain  estimates  of  the  targets' 
locations,  we  use  average,  steady-state,  complex  weight  vectors  in  the  X  and  Y  di¬ 
rections,  that  is. 
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(2.60) 


and 


y 


\r 

m=— M' 

Note  that  each  of  the  P  x  1  average,  steady-state,  complex  weight  vectors  in  the  X  and 
Y  directions  is  given  by 


4r Jl)  =  KJq.l) . KJq.Prf 


(2.62) 


and 


kyjq)  =  fy.  Jq,  1),  hyJq,2). 


(2.63) 


Next,  if  we  now  use  the  average,  steady-state,  complex  weight  vectors  in  the  X  and 
Y  directions,  then  the  one-step,  steady-state,  forward  linear  prediction  errors  for  har¬ 
monic  number  q  in  the  X  and  Y  directions,  are  given  by 


r 

FP'{qJ,n)  =  R{q.i,n)  -  y  ng{q,p)R{qJ  +  p,n),  n  =  -  A'\...,0,...,.Y\  (2.64) 


p=  i 


and 


r 

FPf{q,m,i)  =  R{q.m,i)  -  y  ^hjtg{q,m.p)R{q,m,i  +  p),  m  =  -  M',..., 0,...,Af.  (2.65) 


If  we  define  the  forward  spatial  Z-transform  of  g(i)  as 
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(2.66) 


oc 

G(-)  =  Z{g(i)}  =  £g(/)z', 

i~~oo 


then 


Z{g(i  +  p)}  =  G(z)z~p,  (2.67) 

and,  upon  taking  the  forward  spatial  Z-transform  of  Eqs.  (2.64)  and  (2.65)  with  respect 
to  i,  we  obtain 


(2.68) 


and 


Fp{qjn.zy)  =  R(q,m,zy)[l  - 


m  =  -  M' . 0 . .If. 


(2.69) 


From  Eqs.  (2.68)  and  (2.69),  the  transfer  functions  of  the  one-step  forward  linear 
prediction-error  filters  (analysis  filters)  in  the  X  and  Y  directions  are  given  by 


F i p  (q.zx.n)  I  v—i 

»PE^n)  =  -j—r  =1-4*  2J^'PKP'  n  =  -  A''-° A”’  (2-70> 


R{q,zx.n) 


P=> 


and 


FJq.m}zv) 


p 

HPBy{q.m,Zy)  =  2(q.m,z[)  =  1  “  ~P  m  =  ~  Af . '  (2'71) 


p=\ 


Since  the  average  complex  weights  in  the  X  and  Y  directions  are  independent  of  the  in¬ 
dices  n  and  m,  the  transfer  functions  are  also  independent  of  n  and  m.  that  is. 
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(2.72) 


p 

Hpex{<:Jx)  =  1  ~  -J  - 

p=  i 

and 

p 

H,E,{q.2, )  =  1  -  jr  YfrJl'PP?-  P-73) 

P=  1 


Therefore  the  transfer  functions  of  the  inverse  filters  (synthesis  filters)  in  the  X  and  Y 
directions  are  given  by 


and 


(2.74) 


(2.75) 


In  order  to  obtain  the  spatial-frequency  response  (angular  response)  of  the  filter 
described  by  Eq.  (2.74).  let 


=  exp(  +j2nfxdx)  =  exp( 


+j2nudx 


(2.76) 


Substituting  Eq.(2.76)  into  Eq.(2.74)  yields 


(<7-n) 


1  - 


p 

TV* 


p  exP< 

F=  1 


-jin pud x 


(2.77) 
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Similarly,  if  we  let 


+j2nvdv 

Zy  =  exp(  +j2nfydy)  =  exp( - - )  (2.78) 


in  Eq.  (2. "5),  we  obtain 


1 


—j2npvdy 

/ 


) 


(2.79) 


Equations  (2-77)  and  (2-79)  are  then  evaluated  as  functions  of  the  direction  cosines 
u  and  v,  respectively,  for  a  given  harmonic  q.  The  values  of  u  and  v  that  maximize  Eqs. 
(2-77)  and  (2-79),  respectively,  then  become  the  linear  prediction-error  estimates  u£c(q) 
and  v££(q)  of  ua(q)  and  v0(q)  .  The  estimates  of  d0(q)  and  \p0{q)  at  each  harmonic  q  are 
given  by  [Ref.  3] 


d0(q)  =  sin-'[{[ \uPaE(q)f  +  [tf£(?)]2}'/2],  q  *  0  (2.80) 

and 


•Ao  (q)  =  ^n 


AP£,  .. 
v0  (?) 
'•PE,  \ 
«0  (?) 


q  #  0 


(2.81) 


where  ulE{q)  and  i^E(q)  are  the  estimates  of  the  direction  cosines  obtained  by  using  the 
adaptive  linear  prediction-error  algorithm. 
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III.  SIMULATION  RESULTS 


Computer  simulation  results  for  five  test  cases  are  presented  in  this  chapter.  The  test 
cases  were  designed  to  test  both  of  the  algorithms  multiple  broadband  target  localization 
capability,  full  angular  coverage  capabilities,  and  their  angular  resolution  as  a  function 
of  the  input  SNR  at  a  single  element  in  the  array,  harmonic  number,  and  the  number 
of  iterations. 

For  simplification  purposes,  in  all  of  the  simulation  results  we  use  the  following 
notation: 

•  LMS:  When  the  frequency  domain  modified  complex  Least-Mean-Square  adaptive 
beamforming  algorithm  is  used. 


•  ALP:  When  the  frequency  domain  modified  complex  adaptive  linear  prediction- 
error  filtering  beamforming  algorithm  is  used. 

The  test  cases  are  described  below: 

1.  Case  1  -  A  single  broadband  target  radiating  six  harmonics,  located  at  broadside 
relative  to  the  planar  array; 

2.  Case  2  -  Three  broadband  targets,  each  radiating  two  unique  harmonics,  located 
at  random  positions; 

3.  Case  3  -  A  single  broadband  target  radiating  six  harmonics,  located  at  endfire  rel¬ 
ative  to  the  planar  array; 

4.  Case  4  -  Three  broadband  targets.  Two  targets  are  located  in  the  same  plane  and 
have  a  common  harmonic.  Three  subcases  were  examined.  Each  subcase  used  a 
different  harmonic  as  the  common  spectral  line. 

5.  Case  5  -  Three  broadband  targets,  located  at  random  positions.  Two  targets  have 
a  common  harmonic.  Three  subcases  were  examined.  Each  subcase  used  a  different 
harmonic  as  the  common  spectral  line. 

The  simulation  results  of  these  five  cases  are  based  upon  processing  the  output 
electrical  signals  from  a  1 1  x  1 1  planar  array  of  equally  spaced  hydrophones.  The 
acoustic  field  incident  upon  the  planar  array  was,  in  general,  the  sum  of  several  plane- 
wave  fields  travelling  in  different  directions.  Each  plane-wave  field  consisted  of  an  arbi¬ 
trary  number  of  harmonics  (spectral  lines)  emanating  from  one  of  the  broadband  sound 
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sources.  As  a  result,  the  output  electrical  signal  at  each  element  of  the  planar  array  was 
composed  of  an  arbitrary  number  K  of  harmonics,  all  with  identical  amplitudes  of  unity. 

The  fundamental  frequency  for  all  test  cases  was  chosen  to  be  f0  =  1000  Hz.  Via  Eq. 
(2-6),  the  fundamental  period  at  each  element  of  the  array  was  1  millisecond.  The  sam¬ 
pling  parameter,  S,  was  set  equal  to  10.  The  number  of  time  samples  taken  per  element 
of  the  array  is  given  by 


L  =  SK+  5  (3.1) 

where  L  is  the  total  number  of  samples  taken  per  element,  K  is  the  total  number  of 
harmonics  present  in  the  signal,  and  S  is  the  sampling  parameter.  With  S=  10  and  K  =  6, 
L  =  65  time  samples  per  element. 

In  the  computer  simulations  presented  in  this  thesis,  the  interelement  spacing  is 


where  the  minimum  wavelength  is 


^■min  c 
J  m 


and 


/max  =  */o  Hz 


Substituting  Eqs.  (3.3)  and  (3.4)  into  Eq.  (3.2)  yields 


(3.2) 


(3.3) 


(3.4) 


(3.5) 


For  the  ALP  algorithm,  the  order  of  the  forward  prediction-error  filter  was  10.  The 
reference  element  was  the  first  element  of  each  linear  array  and  the  number  of  linear 
arrays  that  have  been  processed  in  both  directions  was  11.  Initial  simulation  tests  indi- 
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cated  that  the  number  of  iterations  required  for  the  LMS  algorithm  to  approach  a 
steady-state  error  was  100,  while  for  the  ALP  algorithm  1000  iterations  were  required. 

Baseline,  or  "no  noise"  test  case  results  were  the  first  generated  to  ensure  that  both 
algorithms  were  working  properly  in  a  noise-free  environment.  Figure  3  depicts  the 
noise-free,  time  domain,  received  signal  at  element  (1,1)  for  Case  2.  For  each  baseline 
test  case,  bearing  and  depression  angle  estimation  errors  (measured  in  degrees)  were 
obtained  by  running  the  computer  simulation  once  allowing  100  iterations  for  the  LMS 
algorithm  and  1000  iterations  for  the  ALP  algorithm. 

Following  compilation  of  the  baseline  results  for  the  first  three  cases,  identical  test 
cases  were  run  using  additive,  wide-sense  stationary,  zero  mean,  white,  gaussian  noise 
samples  to  corrupt  the  time  samples  of  the  received  signal.  Figure  4  depicts  the  time 
domain,  received  signal  at  element  (1,1)  for  case  2  for  a  SXR  =  -9dB.  For  each  test  case, 
and  for  a  given  input  SNR  at  a  single  element  of  the  array,  average  bearing  and  de¬ 
pression  angle  estimation  errors  were  obtained  by  running  the  computer  simulation  10 
times.  During  each  run.  the  LMS  algorithm  was  allowed  100  iterations  whereas  the  ALP 
algorithm  was  allowed  1000  iterations. 
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A.  CASE  1 


Case  1  placed  a  single  broadband  target  at  broadside  relative  to  the  planar  array 
(i.e.,  d0  =  0°).  This  was  considered  the  simplest  case  for  both  of  the  algorithms  since  it 
is  at  broadside  that  the  far-field  beam  pattern  has  its  narrowest  beamwidth  and,  as  a 
result,  both  algorithms  should  provide  good  angular  resolution  [Ref.  4].  The  general 
plane-wave  field  radiated  by  the  target  consisted  of  six  harmonics. 

Figures  5,  6  and  7  present  the  average  estimation  errors  of  the  depression  angle 
( eBo )  versus  SNR  for  Case  LMS1  at  1000  Hz,  3000  Hz  and  6000  Hz,  respectively.  Results 
for  bearing  angle  estimation  errors  in  the  broadside  case  are  irrelevant  since  the  target 
is  directly  above  the  array,  and  the  bearing  angle  has  no  meaning.  From  these  figures, 
we  can  see  two  trends.  Increasing  the  SNR  results  in  smaller  magnitudes  of  the  average 
estimation  errors  of  the  depression  angle  and  as  the  harmonic  number(q)  increases,  the 
magnitude  of  the  estimation  error  decreases.  These  two  trends  were  expected.  As  we 
increase  the  SNR,  the  noise  components  of  the  received  signals  are  less  dominant.  An 
increase  of  the  harmonic  number  (q)  represents  an  increase  in  the  frequency  and,  as  a 
result,  the  beamwidth  of  the  far-field  beam  pattern  decreases.  The  decreased  beamwidth 
increases  the  angular  resolution  of  the  algorithm  and,  as  result,  decreases  the  estimation 
error. 

Figures  8,  9  and  10  present  estimates  of  the  direction  cosines  lT0  and  Va  for  Case 
ALP1  at  SNR=-3  dB  and  frequencies  1000  Hz,  3000  Hz,  and  6000  Hz,  respectively. 
Similarly,  Figures  11,  12  and  13,  and  14,  15,  and  16  present  estimates  of  the  direction 
cosines  i'0  and  V0  for  Case  ALP1  with  the  same  order  of  frequency  as  above  and  for 
SN'R  =  -6  dB  and  SNR  =  -9  dB,  respectively.  As  we  discussed  in  Chapter  2,  the  location 
of  the  peaks  of  the  magnitude  of  the  synthesis  filter  transfer  function  in  dB  give  the  U0 
and  J’0  estimates.  From  these  figures  we  can  clearly  observe  that  as  the  SNR  decreases, 
the  difference  between  the  peak  and  the  background  noise  in  the  transfer  function  plots 
decreases.  Also,  the  harmonic  number  is  important,  because  increasing  the  harmonic 
number  produces  sharper  peaks  in  the  transfer  function  plots  resulting  in  more  accurate 
estimation. 

Figures  17,  18  and  19  present  the  average  estimation  error  of  the  depression  angle 
( e9o )  versus  SNR  for  Case  ALP1  at  1000  Hz,  3000  Hz  and  6000  Hz.  The  trends  which 
were  present  in  Case  L.MS1  are  apparent  here,  that  is,  as  SNR  increases,  the  magnitude 
of  the  estimation  error  decreases  and,  as  the  harmonic  number  (q)  increases,  the  mag¬ 
nitude  of  the  estimation  error  decreases. 
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Finally  Tables  1  and  2  in  the  Appendix  A  present  the  numerical  data  for  Cases 
LMS1  and  ALP1  respectively  for  all  six  harmonics.  Comparing  the  two  algorithms  us¬ 
ing  these  tables,  we  can  conclude  that  for  Case  1,  both  algorithms  give  small  magnitude 
estimation  errors  in  low  SNR  with  a  little  better  performance  for  the  LMS1  particularly 
at  1000  Hz. 
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FREQ:  1000.0  HZ  M-  11  N-  11  NITER-  100 


Figure  5.  Average  Estimation  Error  of  Depression  Angle  vs.  SNR  at  1000  Hz  for 
Case  LMS1. 
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Figure  6.  Average  Estimation  Error  of  Depression  Angle  vs.  SNR  at  3000  Hz  for 
Case  LMSl. 
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7.  Average  Estimation  Error  of  Depression  Angle  vs.  SNR  at  6000  Hz  for 
Case  LMSI. 
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Figure  9.  Estimates  of  Direction  Cosines  UO  and  VO,  Case  ALP1:  SNR  =  -3dB  and 
q  =  3 
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Figure  14.  Estimates  of  Direction  Cosines  UO  and  VO,  Case  ALP1:  SNR  =  -9dB  and 
q=  1 
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Figure  17.  Average  Estimation  Error  of  Depression  Angle  vs.  SNR  at  1000  Hz  for 
Case  ALP1. 
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Figure  18.  Average  Estimation  Error  of  Depression  Angle  vs.  SNR  at  3000  Hz  for 
Case  ALP1. 
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Figure  19.  Average  Estimation  Error  of  Depression  Angle  vs.  SNR  at  6000  Hz  for 
Case  ALP1. 
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B.  CASE  2 


Case  2  placed  three  broadband  targets  at  random  locations.  The  general  plane-wave 
field  radiated  by  each  of  the  targets  contained  two  unique  harmonics  .  Target  1  was  lo¬ 
cated  at  (0O  =  49°,  =  38°)  and  radiated  harmonics  1  and  6.  Target  2  was  located  at 

(da  =  5°,  <p0  =  137°)  and  radiated  harmonics  2  and  5.  Target  3  was  located  at  (0O  =  77°, 
i/'o  =  307°)  and  radiated  harmonics  3  and  4.  Since  there  are  three  incident  plane-wave 
fields,  each  containing  two  unique  harmonics,  the  output  electrical  signal  from  each  el¬ 
ement  in  the  array  exhibits  a  total  of  six  harmonics. 

Figures  20  through  25  present  the  average  estimation  errors  of  the  depression  (eSg) 
and  bearing  (e*0)  angles  versus  SNR  for  Case  LMS2  for  all  the  harmonics  in  increasing 
order  of  the  harmonic  number  (q),  respectively.  The  same  trend  about  the  SNR  as  in 
Case  1  holds  here,  that  is  increasing  the  SNR  results  in  smaller  magnitude  of  the  average 
estimation  error.  These  general  results  can  be  explained  using  the  same  arguments  pre¬ 
sented  for  Case  1. 

Figures  26  through  43  present  estimates  of  the  direction  cosines  U0  and  V0  for  all  the 
harmonics  placed  in  increasing  order  of  the  harmonic  number  (q)  at  SNR's  -3  dB,  -6  dB 
and  -9  dB.  respectively.  Again,  the  two  trends  observed  in  Case  ALP1  are  apparent  here, 
that  is.  decreasing  the  SNR  decreases  the  difference  between  the  peak  and  the  back¬ 
ground  noise  and  increasing  the  harmonic  number  (q)  produces  sharper  peaks  in  the 
transfer  function  plots. 

Figures  44  through  49  present  the  average  estimation  errors  of  the  depression  (eBn) 
and  bearing  (t\0)  angles  versus  SNR  for  Case  ALP2  for  all  the  harmonics  in  increasing 
order  of  the  harmonic  number  (q).  respectively.  The  same  trends  as  in  LMS2  hold,  with 
the  exception  that  now  the  magnitude  of  the  average  estimation  error  for  both  angles  is 
smaller  than  in  LMS2  .  Note  the  difference  in  scaling  between  the  error  plots  for  Case 
LMS2  and  Case  ALP2. 

Finally,  Tables  3  and  4  in  the  Appendix  A  present  the  numerical  data  for  Cases 
LMS2  and  ALP2,  respectively,  for  all  six  harmonics.  Comparing  the  two  algorithms 
using  these  tables,  we  can  conclude  that  for  Case  2,  the  ALP  algorithm  gives  smaller 
magnitude  estimation  errors  in  most  instances  than  the  LMS  algorithm  especially  for 
low  SNR. 
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Figure  20.  Average  Estimation  Errors  of  Depression  and  Bearing  Angles  vs.  SNR 
at  1000  Hz  for  Case  LMS2. 
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Figure  21.  Average  Estimation  Errors  of  Depression  and  Bearing  Angles  vs. 
at  2000  Hz  for  Case  LMS2. 
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Figure  22.  Average  Estimation  Errors  of  Depression  and  Bearing  Angles  vs.  SNR 
at  3000  Hz  for  Case  LMS2. 
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Figure  23.  Average  Estimation  Errors  of  Depression  and  Bearing  Angles  vs.  SNR 
at  <1000  Hz  for  Case  LMS2. 


47 


AVERAGE  ESTIMATION  ERROR (DEG) 

-30.0  -24.0  -18.0  -12.0  -6.0  0.0  6.0  12.0  1B.0  24.0  30.0 


-9.0  -6.0  -3.0 


input  snr (DO)  rrr  n  single  drrdy  elemcnt 

CASE: LMS2 

FREO:  5000.0  HZ  M-  1 1  N-  1 1  NITER-  100 

Figure  24.  Average  Estimation  Errors  of  Depression  and  Bearing  Angles  vs.  SNR 
at  5000  Hz  for  Case  LMS2. 
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Figure  25.  Average  Estimation  Errors  of  Depression  and  Bearing  Angies  vs.  SNR 
at  6000  Hz  for  Case  LMS2. 
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Figure  26.  Estimates  of  Direction  Cosines  UO  and  VO,  Case  ALP2:  SNR  =  -3dB  and 
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Figure  28.  Estimates  of  Direction  Cosines  Ul)  and  YU,  Case  ALP2:  SNR  =  -3dB  and 
q=3 
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Figure  33.  Estimates  of  Direction  Cosines  UO  and  VO,  Case  ALP2:  SNR  =  -6dB  and 
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Figure  35.  Estimates  of  Direction  Cosines  UU  and  VU,  Case  ALF2:  SNR  =  -6dB  and 
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Figure  36.  Estimates  of  Direction  Cosines  UO  and  VO,  Case  ALP2:  SNR  =  -6dB  and 
q  =  5 
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Figure  38.  Estimates  of  Direction  Cosines  Lit)  and  VO,  Case  ALP2:  SNR  =  -9dB  and 
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Figure  41.  Estimates  of  Direction  Cosines  UO  and  VO,  Case  ALP2:  SNR  =  -9dB  and 
q  =  4 
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Figure  42.  Estimates  of  Direction  Cosines  UO  and  VO,  Case  ALP2:  SNR  =  - 
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Figure  43.  Estimates  of  Direction  Cosines  UU  and  VO,  Case  ALP2:  SNR  =  -9dB  and 
q  =  6 
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AVERAGE  ESTIMATION 

-10.0  -8.0  -6.0  -4.0  -2.0  0.0  2.0  4 


-9.0  -6.0  -3.0 


INPUT  SNR(DQ)  AT  A  SINGLE  ARRAY  ELEMENT 

CASE: ALP2 

EREQ:  1000.0  HZ  M-  1 1  N-  1 1  I-  1  P-  10  NITER-  1000 


Figure  44.  Average  Estimation  Errors  of  Depression  and  Bearing  Angles  vs.  SNR 
at  1000  Hz  for  Case  ALP2. 
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AVERAGE  ESTIMATION  ERROR(OEG) 
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Figure  45.  Average  Estimation  Errors  of  Depression  and  Bearing  Angles  vs.  SNR 
at  2000  Hz  for  Case  ALP2. 
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RVERRGE  ESTIMATION  ERROR  (DEG) 
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Figure  46.  Average  Estimation  Errors  of  Depression  and  Bearing  Angles  vs.  SNR 
at  3000  Hz  for  Case  ALP2. 
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RVERRGE  ESTIMATION  ERROR (DEG) 
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Figure  47.  Average  Estimation  Errors  of  Depression  and  Bearing  Angles  vs.  SNR 
at  4000  Hz  for  Case  ALP2. 
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Figure  49.  Average  Estimation  Errors  of  Depression  and  Bearing  Angles  vs.  SNR 
at  6000  Hz  for  Case  ALP2. 
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C.  CASE  3 


Case  3  placed  a  single  broadband  target  at  endfire  relative  to  the  planar  array  (i.e.. 
60  =  90°,  t//0  =  90°).  This  case  was  considered  the  most  difficult  for  both  algorithms  since 
it  is  at  endfire  that  the  far-field  beam  pattern  beamwidth  is  broadest  [Ref.  4].  The  general 
plane-wave  field  radiated  by  the  target  consisted  of  six  harmonics. 

Figures  50  and  51  present  the  average  estimation  errors  of  the  depression  (e<,0)  and 
bearing  (e„0)  angles  versus  SNR  for  Case  LMS3  at  1000  Hz  and  3000  Hz,  respectively. 
The  estimation  errors  for  the  highest  harmonic  (6000  Hz)  are  not  shown  for  Case  LMS3 
since  the  magnitude  of  these  errors  were  between  45  and  90  degrees.  This  poor  per¬ 
formance  for  the  highest  harmonic  at  endfire  can  be  explained  theoretically  [Ref.  6]. 

Figures  52,  53,  and  54  present  estimates  of  the  direction  cosines  UD  and  F0  for  Case 
ALP3  at  SNR=-3  dB  and  frequencies  1000  Hz,  3000  Hz,  and  6000  Hz,  respectively. 
Similarly,  Figures  55,  56,  and  57,  and  58.  59,  and  60  present  estimates  of  the  direction 
cosines  L'0  and  V0  for  Case  ALP3  with  the  same  order  of  frequency  as  above  and 
SNR=-6  dB  and  SNR=-9  dB.  respectively.  Here  we  can  observe  an  ambiguity  for  the 
F0  estimates  at  6000  Hz  which  can  be  explained  by  examining  the  exponential  term  of 
Eq.(2-79).  This  term  for  the  highest  frequency  component  creates  integer  multiplies  of 
2n  producing  two  peaks  at  ±  1  for  the  V0  estimates  at  6000  Hz. 

Figures  61,  62  and  63  present  the  average  estimation  errors  of  the  depression  (ce ) 
and  bearing  (evJ  angles  versus  SNR  for  Case  ALP3  at  1000  Hz,  3000  Hz,  and  6000  Hz, 
respectively.  We  can  observe  that  in  Case  ALP3  at  6000  Hz  we  have  acceptable  esti¬ 
mates  for  depression  angle  (6C).  In  Figure  63  the  \p0  error  plot  is  not  shown  since  the 
magnitude  of  this  error  was  between  45  and  90  degrees. 

Finally,  Tables  5  and  6  in  Appendix  A  present  the  numerical  data  for  Cases  LV1S3 
and  ALP3  respectively  for  all  six  harmonics.  Comparing  the  two  algorithms  using  these 
tables,  we  can  conclude  that  for  Case  3  the  ALP  algorithm  gives  smaller  magnitude  es¬ 
timation  errors  in  most  instances  than  the  LMS  algorithm  for  both  angles. 
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Figur;  50.  Average  Estimation  Errors  of  Depression  and  Bearing  Angies  vs.  SNR 
at  1000  Hz  for  Case  LMSJ. 
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Figuio  51.  Average  Estimation  Errors  of  Depression  and  Bearing  Angles  vs.  SNR 
at  3000  Hz  for  Case  LMS3. 
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Figure  54. 


CflSE:flLP3  SI*-  -3.0  00 

FREOi  6000.0  HZ  M-  It  N-  1 1  I-  I  P-  10  NITER-  1000 
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Figure  55.  Estimates  of  Direction  Cosines  UO  and  VO,  Case  ALP3:  SNR  =  -6dB  and 
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Figure  56.  Estimates  of  Direction  Cosines  UU  and  VO,  Case  ALF3:  SNR  =  -6dB  and 
q  =  3 
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Figure  57.  Estimates  of  Direction  Cosines  UO  and  VO,  Case  ALP3:  SNR  =  -6dB  and 
q  =  6 
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CASE: BLP3 

FREQ:  1000.0  HZ  M-  1 1  N-  1 1  I-  1  P-  10  NITER-  1000 

Figure  6i.  Average  Estimation  Errors  of  Depression  and  Bearing  Angles  vs.  SNR 

at  1000  Hz  for  Case  ALP3. 
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Figure  62.  Average  Estimation  Errors  of  Depression  and  Bearing  Angies  vs.  SNR 
at  3000  Hz  for  Case  ALP3. 
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Figure  63.  Average  Estimation  Errors  of  Depression  and  Bearing  Angles  vs.  SNR 
at  6000  Hz  for  Case  ALP3. 


88 


D.  CASE  4 


Case  4  corresponds  to  three  broadband  targets  being  present  two  of  which  have  a 
common  harmonic.  Target  1  was  located  at  ( 60  =  45°,  =  0°).  Target  2  was  located  at 

(0O  =  45°,  =  180°).  Target  3  was  located  at  (  60  =  33°  ,  =  47°).  A  total  of  six  har¬ 

monics  were  present  in  the  output  electrical  signal.  The  following  three  subcases  were 
examined: 

•  Case  4-A:  Three  targets,  two  of  which  are  located  in  the  same  plane.  Targets  1  and 
2  have  a  common  spectral  line  at  1000  Hz. 

•  Case  4-B:  Three  targets,  two  of  which  are  located  in  the  same  plane.  Targets  1  and 
2  have  a  common  spectral  line  at  3000  Hz. 

•  Case  4-C:  Three  targets  two  of  which  are  located  in  the  same  plane.  Targets  1  and 
2  have  a  common  spectral  line  at  6000  Hz. 

For  these  cases,  only  no  noise  results  were  compiled  for  both  algorithms  in  order  to 
examine  the  ability  of  target  discrimination  when  two  targets  are  in  the  same  plane  and 
have  a  common  harmonic. 

Figures  64,  65  and  66  present  estimates  of  the  direction  cosines  U0  and  VB  for  each 
common  harmonic,  that  is  at  1000  Hz  for  Case  ALP4-A,  3000  Hz  for  Case  ALP4-B  and 
at  6000  Hz  for  Case  ALP4-C.  Clearly,  we  can  observe  that  except  for  case  ALP4-A,  the 
algorithm  can  discriminate  the  location  of  two  targets. 

Tables  7  through  12  in  the  Appendix  A  present  numerical  data  for  all  the  subcases 
of  Case  4  for  both  algorithms.  From  these  tables  we  can  observe  that  for  all  the  har¬ 
monics,  except  the  common  one.  both  algorithms  correctly  identified  the  location  of  the 
targets  with  zero  degrees  estimation  errror.  However,  for  the  common  harmonic,  the 
LMS  algorithm  located  a  false  fourth  target,  which  was  located  exactly  between  the  ac¬ 
tual  locations  of  targets  1  and  2.  On  the  other  hand  the  ALP  algorithm,  except  for  the 
subcase  ALP4-A,  for  the  common  harmonic,  identifies  two  different  estimates  of  the  di¬ 
rection  cosine  C'0  from  the  transfer  function  plots,  which  results  in  two  different  target 
locations. 
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Figure  64.  Estimates  of  Direction  Cosines  UO  and  VO,  Case  ALP4-A:  No  Noise  and 

q=  i 
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Figure  65.  Estimates  of  Direction  Cosines  UO  and  VO,  Case  ALP4-B:  No  Noise  and 
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Figure  66.  Estimates  of  Direction  Cosines  UO  and  VO,  Case  ALP4-C:  No  Noise  and 


E.  CASE  5 


Case  5  placed  three  broadband  targets  at  random  locations.  Target  1  was  located 
at  ( 80  =  49°,  i l/a  =  38°).  Target  2  was  located  at  ( 60  =  5°,  \p0  =  137°).  Target  3  was  located 
at  <0O  =  77°,  \p0  =  307°).  A  total  of  six  harmonics  were  present  in  the  output  electrical 
signal.  The  following  three  subcases  were  examined: 

•  Case  5-A:  Three  targets  located  at  random  positions.  Targets  I  and  2  have  a 
common  spectral  line  at  1000  Hz. 

•  Case  5-B:  Three  targets  located  at  random  positions.  Targets  1  and  2  share  a 
common  spectral  line  at  3000  Hz. 

•  Case  5-C :  Three  targets  located  at  random  positions.  Targets  1  and  2  share  a 
common  spectral  line  at  6000  Hz. 

For  these  cases,  only  no  noise  results  were  compiled  for  both  algorithms  in  order  to 
examine  the  ability  of  target  discrimination  when  two  targets  are  in  random  positions 
and  have  a  common  harmonic. 

Figures  67.  68  and  69  present  estimates  of  the  direction  cosines  i'0  and  V0  for  each 
common  harmonic,  that  is.  at  1000  Hz  for  Case  ALP5-A,  at  3000  Hz  for  Case  ALP5-B 
and  at  6000  Hz  for  case  ALP5-C.  Clearly,  we  can  observe  that  exept  for  case  ALP5-A, 
the  algorithm  can  discriminate  the  location  of  two  targets. 

Tables  13  through  18  in  the  Appendix  A  present  numerical  data  for  all  the  subcases 
of  Case  5  for  both  algorithms.  From  these  tables  we  can  observe  that  for  all  the  har¬ 
monics,  except  the  common  one,  both  algorithms  correctly  identified  the  location  of  the 
targets  with  zero  degrees  estimation  errror.  However,  for  the  common  harmonic,  the 
LMS  algorithm  located  a  false  fourth  target.  On  the  other  hand  the  ALP  algorithm, 
except  for  subcase  ALP5-A.  for  the  common  harmonic,  identifies  two  different  estimates 
of  the  direction  cosines  U0  and  V0  from  the  transfer  function  plots,  which  results  in  two 
different  target  locations. 
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Figure  67.  Estimates  of  Direction  Cosines  UO  and  VO,  Case  ALP5-A:  No  Noise  and 
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Figure  69.  Estimates  of  Direction  Cosines  UU  and  VO,  Case  ALP5-C:  No  Noise  and 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  goal  of  this  thesis  was  to  evaluate  the  multiple  broadband  target  localization 
capability  and  the  full  angular  coverage  capability  of  two  frequency  domain  adaptive 
beamforming  algorithms.  The  main  advantage  of  both  algorithms  is  that  no  a  priori 
knowledge  is  assumed  about  the  statistics  of  the  received  signal,  the  total  number  of 
targets  present,  where  they  are  located,  and  the  frequency  content  of  the  targets. 

The  two  frequency  domain  adaptive  beamforming  algorithms  evaluated  were  the 
modified  complex  Least-Mean-Square  (LMS)  algorithm  developed  by  Ziomek  and  Chan 
[Ref.  6],  and  the  modified  complex  Adaptive  Linear  Prediction-Error  (ALP)  algorithm 
which  was  presented  for  first  time.  Chapter  3  presented  the  results  of  five  test  cases 
which  were  designed  to  test  both  algorithms  for  several  different  capabilities.  It  should 
also  be  noted  here  that  the  number  of  time  samples  taken  per  element  of  the  array  never 
exceeded  65  time  samples,  an  amount  which  represents  a  small  number  of  data  points. 

Summarizing  the  simulation  results  discussed  in  Chapter  3,  we  can  list  the  following 
conclusions'. 

•  In  most  cases  tested,  both  algorithms  can  localize  multiple  broadband  targets  and 
they  have  full  angular  coverage  capability. 

•  Increasing  either  the  sampling  rate  or  the  array  size  yields  better  estimates  in  lower 
SNR  for  both  algorithms 

•  ALP  generally  gives  better  estimates  in  lower  SNR  than  LMS 

•  ALP  can,  in  most  of  the  "no  noise"  cases,  discriminate  between  two  targets  that 
have  a  common  harmonic  while  LMS  cannot 

•  ALP  is  more  computationally  intensive  than  LMS  and  requires  much  more  iter¬ 
ations  to  approach  a  steady-state  error. 

•  ALP  is  more  flexible  than  LMS  because  it  can  process  only  a  part  of  the  total 
number  of  array  elements  to  localize  the  target 

In  the  course  of  this  investigation,  several  possible  areas  for  future  research  pre¬ 
sented  themselves: 
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•  further  study  of  the  effects  of  varying  the  step-size  parameter,  n  ,  in  an  attempt  to 
decrease  the  magnitude  of  the  estimation  error. 

•  investigation  of  parallel  processing  implementation  for  the  ALP  algorithm  in  order 
to  reduce  the  computation  time. 

•  investigation  of  larger  array  size  in  lower  SNR  mainly  using  the  ALP  algorithm 

•  investigation  of  a  noise  reduction  system  prior  to  processing  the  output  electrical 
signals. 


% 
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APPENDIX  A. 

This  appendix  presents  the  numerical  data  for  all  of  the  test  cases.  Tables  1  through 
6  present  numerical  data  for  the  three  first  "noisy"  cases,  Cases  1  through  3,  for  both 
algorithms,  while  tables  7  through  12  present  numerical  data  from  the  two  "no  noise" 
cases,  Case  4  and  5,  for  both  algorithms.  The  following  notation  is  used  in  all  the  tables: 

•  LMS:  When  the  frequency  domain  modified  complex  Least-Mean-Square  adaptive 
beamforming  algorithm  is  used. 

•  ALP:  When  the  frequency  domain  modified  complex  adaptive  linear  prediction- 
error  filter  beamforming  algorithm  is  used 

•  60  :  Depression  Angle 

•  i^0  :  Bearing  Angle 

•  c„o :  Average  Estimation  Error  for  Depression  Angle 

•  evg  :  Average  Estimation  Error  for  Bearing  Angle 
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Table  1.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  LMS1 


Case  1:  One  target  located  at  broadside  relative  to  planar  array  -  six  harmonics 
present  in  output  electrical  signal  -  LMS  algorithm  is  used  -  100  iterations 


Average  Estimation  Error  (deg) 


SNR 


-6  dB 


Harmonic  (q) 


Location  (deg) 

lAo 

0 

270 

0 

270 

0 

270 

0 

270 

0 

270 

0 

270 

-3  dB 


e 


-2.169 


-0.810 


-0.561 


-0.360 


-0.362 


-0.299 


-1.789 


-1.426 


-0.607 


-0.760 


-0.363 


-0.435 


-9  dB 


e 


-4.382 


-2.170 


-1.233 


-0.649 


-0.S48 


-1.596 


Table  2.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  ALP1 


Case  1:  One  target  located  at  broadside  relative  to  planar  array  -  six  harmonics 
present  in  output  electrical  signal  -  ALP  algorithm  is  used  -  1000  iterations 


Harmonic  (q) 


Location  (deg) 

0o 

<Ao 

0 

270 

0 

270 

0 

270 

0 

270 

0 

270 

0 

270 

Average  Estimation  Error  (deg) 


SNR 


-6  dB 


-3  dB 


-3.623 


-1.307 


-1.325 


-1.231 


-0.753 


-0.471 


-4.444 


-2.739 


-1.831 


-1.429 


-0.740 


-0.665 


-9  dB 


*»« 


-6.056 


-2.404 


-1.993 


-1.463 


-1.746 


-0.951 


Table  3.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  LMS2 


Case  2:  Three  targets  located  at  random  positions  -  six  harmonics  present  in  output 
electrical  signal  -  LMS  algorithm  is  used  -  100  iterations 


Harmonic 

(q) 

Location  (deg) 

Average  Estimation  Error  (deg) 

e„ 

to 

SNR 

-3  dB 

-6  dB 

-9  dB 

e°n 

An 

e°a 

e*n 

^n 

<V 

i 

49 

38 

0.265 

-0.124 

-0.411 

-0.832 

-0.571 

0.325 

2 

5 

137 

0.077 

-2.714 

-0.161 

-5.017 

0.097 

3.451 

■> 

-3 

77 

307 

3.037 

-0.108 

0.351 

-2.625 

7.991 

9.333 

4 

77 

307 

-10.532 

1.677 

-5.032 

0.019 

3.858 

26.436 

5 

5 

137 

-0.022 

0.560 

-0.155 

1.376 

0.07S 

0.002 

6 

49 

38 

1.081 

-1.531 

-10.852 

4.382 

-2.355 

-8.784 

Table  4.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  ALP2 


Case  2:  Three  targets  located  at  random  positions  -  six  harmonics  present  in  output 
electrical  signal  -  ALP  algorithm  is  used  -  1000  iterations _ 


Harmonic 

(q) 

Location  (deg) 

Average  Estimation  Error  (deg) 

I 

1 

to 

SNR 

-3  dB 

-6  dB 

-9  dB 

WBM 

_  .  e*n 

0>n 

1 

49 

38 

-1.470 

0.098 

-0.826 

0.380 

8.243 

-2.443 

2 

5 

137 

0.134 

6.094 

-1.091 

-8.072 

0.879 

1.252 

3 

77 

307 

-1.72S 

0.082 

-3.247 

-0.361 

-0.697 

1.131 

4 

77 

307 

-1.968 

0.526 

-0.964 

0.21S 

-1.763 

-0.035 

5 

5 

137 

0.019 

-0.191 

0.180 

1.267 

0.830 

-6.591 

6 

49 

38 

-0.173 

0.146 

-0.197 

-0.019 

2.06S 

-5.109 
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Table  5.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  LMS3 


Case  2:  One  target  located  at  endfire  relative  to  planar  array  -  six  harmonics 
in  output  electrical  signal  -  LMS  algorithm  is  used  -  100  iterations 


Location  (deg) _ Average  Estimation  Error  (deg) 

Harmonic  _ _ SNR _ 

(q)  en  I h  -3  dB  -6  dB 


present 


Location  (deg) 

00 

'/'o 

90 

90 

90 

90 

90 

90 

90 

90 

90 

90 

90 

90 

3.355  -0.055 


33.033  -0.200 


19.490  -0.019 


2.247  0.006 


2.136  -0.0S3 


70.S77  -84.120  69.426 


20.000 


40.425 


16.401 


7.365 


1.146 


0.545 


-0.646 


-0.229 


0.287 


0.09S 


-64.811 


-9  dB 


e 


21.186 


27.226 


16.706 


12.025 


14.692 


46.737 


1.023 


-0.330 


-0.325 


-0.139 


-18.528 


-90.393 


Table  6.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  ALP3 


Case  2:  One  target  located  at  endfire  relative  to  planar  array  -  six  harmonics  present 
in  output  electrical  signal  -  ALP  algorithm  is  used  -  1000  iterations 


Harmonic 

(q) 

Average  Estimation  Error  (deg) 

'h 

SNR 

-3  dB 

-6  dB 

-9  dB 

% 

^'n 

*«» 

_  evr\ 

i 

90 

90 

3.948 

0.217 

7.806 

-0.136 

6.672 

-0.984 

2 

90 

90 

7.287 

-0.362 

11.323 

-0.48S 

9.034 

0.329 

3 

90 

90 

5.342 

-0.107 

-0.370 

1.935 

-0.767 

4 

90 

90 

2.233 

-0.113 

2.145 

0.286 

7.725 

-0.023 

5 

90 

90 

0.986 

-0.160 

3.512 

0.076 

14.090 

-6.512 

6 

90 

90 

5.229 

-89.997 

7.051 

-90.029 

16.548 

-S9.0S8 
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Tahle  7.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  LMS4-A 


Case  4a:  Three  targets  two  of  which  are  located  in  the  same  ^'ane  -  targets  one  and 
two  have  a  common  spectral  line  at  1000  Hz  •  six  harmonk  .  piv..ent  in  output  elec¬ 
trical  signal  -  LMS  algorithm  is  used  -  100  iterations  -  no  noise 

Harmonic  (q) 

Location  (deg) 

Average  Estimation  Error 
(deg) 

e0 

1 

45 

0 

45.000 

-90.000 

1 

45 

180 

45.000 

90.000 

2 

45 

0 

0.000 

0.000 

3 

33 

47 

0.000 

0.000 

4 

33  | 

47 

0.000 

0.000 

5 

45 

180 

0.000 

0.000 

6 

33 

47 

0.000 

0.000 

Table  8.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  ALP4-A 


Case  4a:  Three  targets  two  of  which  are  located  in  the  same  plane  -  targets  one  and 
two  have  a  common  spectral  line  at  1000  Hz  -  six  harmonics  present  in  output  elec¬ 
trical  signal  -  ALP  algorithm  is  used  -  1000  iterations  -  no  noise 

Harmonic  (q) 

1 

Location  (deg) 

Average  Estimation  Error 
(deg) 

00 

<Ao 

** 

*Vn 

1 

45 

0 

-44.921 

0.0 

1 

45 

ISO 

-44.921 

ISO.OOO 

*7 

i- 

45 

0 

0.009 

0.000 

3 

33 

47 

0.037 

-0.011 

4 

33 

47 

0.037 

-0.01 1 

5 

45 

180 

0.009 

0.000 

6 

33 

47 

0.037 

-0.011 
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Table  9,  NUMERICAL  DATA  CORRESPONDING  TO  CASE  LMS4-B _ 

Case  4b:  Three  targets  two  of  which  are  located  in  the  same  plane  -  targets  one  and 
two  have  a  common  spectral  line  at  3000  Hz  •  six  harmonics  present  in  output  elec- 
trical  signal  -  LMS  algorithm  is  used  -  100  iterations  -  no  noise _ 

Location  (deg)  |  Average  Estimation  Error 


Harmonic  (q) 


e0 

33 

47 

0.000 

0.000 

45 

0 

0.000 

0.000 

45 

0 

45.00 

-90.00 

45 

180 

45.0 

90.0 

33 

47 

0.000 

0.000 

45 

180 

0.000 

0.000 

33 

47 

0.000 

0.000 

Table  10.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  ALP4-B 


Case  4b:  Three  targets  two  of  which  are  located  in  the  same  plane  -  targets  one  and 
two  have  a  common  spectral  line  at  3000  Hz  -  six  harmonics  present  in  output  elec¬ 
trical  signal  -  ALP  algorithm  is  used  -  1000  iterations  -  no  noise 


Estimation  Error  (deg) 


e 


Location  (deg) 

e0 

<Pn 

33 

47 

45 

0 

45 

0 

45 

180 

33 

47 

45 

ISO 

33 

47 

0.037 


0.009 


1.56 


1.56 


0.032 


0.009 


0.037 


-0.011 


0.000 


0.000 


0.000 


-0.011 


0.000 


-0.011 


Table  11.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  LMS4-C 


Case  4c:  Three  targets  two  of  which  are  located  in  the  same  plane  -  targets  one  and 
two  have  a  common  spectral  line  at  6000  Hz  -  six  harmonics  present  in  output  elec¬ 
trical  signal  -  LMS  algorithm  is  used  -  100  iterations  -  no  noise 

Harmonic  (q) 

Location  (deg) 

Average  Estimation  Error 
(deg) 

ea 

iAo 

1 

33 

47 

0.000 

0.000 

2 

45 

0 

0.000 

0.000 

3 

33 

47 

0.000 

0.000 

4 

33 

47 

0.000 

0.000 

5 

45 

180 

0.000 

o.ooo 

6 

45 

0 

45.000 

-90.000 

6 

45 

ISO 

45.000 

90.000 

Table  12.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  ALP4-C 


Case  4c:  Three  targets  two  of  which  are  located  in  the  same  plane  -  targets  one  and 
two  have  a  common  spectral  line  at  6000  Hz  -  six  harmonics  present  in  output  elec¬ 
trical  signal  -  ALP  algorithm  is  used  -  1000  iterations  -  no  noise 

Harmonic  (q) 

Location  (deg) 

Average  Estimation  Error 
(deg) 

eB 

•Ao 

e»» 

- 

i 

33 

47 

0.037 

-0.011 

2 

45 

0 

0.009 

0.000 

3 

33 

47 

0.037 

-0.011 

4 

33 

47 

0.037 

-0.011 

5 

45 

180 

0.009 

0.000 

6 

45 

0 

2.070 

0.000 

6 

45 

ISO 

2.070 

0.000 
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Table  13.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  LMS5-A 


Case  5a:  Three  targets  located  at  random  positions  -  targets  one  and  two  have  a 
common  spectral  line  at  1000  Hz  -  six  harmonics  present  in  output  electrical  signal  - 
LMS  algorithm  is  used  -  100  iterations  -  no  noise 

Harmonic  (q) 

Location  (deg) 

Average  Estimation  Error 
(deg) 

e„ 

V'o 

e0(t 

e*n 

1 

49 

38 

-13.000 

252.594 

1 

5 

137 

-41.000 

-16.406 

2 

5 

137 

0.000 

0.000 

3 

77 

307 

0.000 

0.000 

4 

77 

307 

0.000 

0.000 

5 

5 

137 

0.000 

0.000 

6 

49 

38 

0.000 

0.000 

Table  14.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  ALP5-A 


Case  5a:  Three  targets  located  at  random  positions  -  targets  one  and  two  have  a 
common  spectral  line  at  1000  Hz  -  six  harmonics  present  in  output  electrical  signal  - 
ALP  algorithm  is  used  -  1000  iterations  -  no  noise 

Harmonic  (q) 

Location  (deg) 

Average  Estimation  Error 
(deg) 

lAo 

C-o 

1 

49 

38 

27.121 

-6.674 

1 

5 

137 

-16.879 

92.326 

2 

5 

137 

-0.038 

-0.008 

3 

77 

307 

0.094 

0.012 

4 

77 

307 

0.094 

0.012 

5 

5 

137 

0.006 

-0.328 

6 

49 

38 

0.094 

0.012 
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Table  15.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  LMS5-B 


Case  5b:  Three  targets  located  at  random  positions  -  targets  one  and  two  have  a 
common  spectral  line  at  3000  Hz  -  six  harmonics  present  in  output  electrical  signal  - 
LMS  algorithm  is  used  -  100  iterations  -  no  noise 

Harmonic  (q) 

Location  (deg) 

Average  Estimation  Error 
(deg) 

0 n 

<Ao 

«•« 

1 

77 

307 

0.000 

0.000 

2 

5 

137 

0.000 

0.000 

3 

49 

38 

-13.000 

252.594 

3 

5 

137 

-41.000 

-16.406 

4 

77 

307 

0.000 

0.000 

5 

5 

137 

0.000 

0.000 

6 

49 

38 

0.000 

0.000 

Table  16.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  ALP5-B 


Case  5b:  Three  targets  located  at  random  positions  - 

targets  one  and  two  have  a  common  spectral  line  at  3000  Hz  -  six  harmonics  present 
in  output  electrical  signal  -  ALP  algorithm  is  used  -  1000  iterations  -  no  noise 

Harmonic  (q) 

Location  (deg) 

Estimation  Error  (deg) 

^'0 

1 

77 

307 

0.094 

0.012 

2 

5 

137 

-0.038 

-0.008 

49 

38 

2.94 

0.39 

3 

5 

137 

1.505 

-1.992 

4 

77 

307 

0.094 

0.012 

5 

5 

137 

0.006 

-0.328 

6 

49 

38 

0.094 

0.012 
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Table  17.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  LMS5-C 


Case  5c:  Three  targets  located  at  random  positions  -  targets  one  and  two  have  a 
common  spectral  line  at  6000  Hz  -  six  harmonics  present  in  output  electrical  signal  - 
LMS  algorithm  is  used  -  100  iterations  -  no  noise 


Average  Estimation  Error 
(deg) 


Harmonic  (q) 


Location  (deg) 

0o 

77 

307 

5 

137 

49 

38 

77 

307 

5 

137 

5 

137 

49 

38 

e 


0.000 


0.000 


0.000 


0.000 


0.000 


-41.000 


-13.000 


0.000 


0.000 


0.000 


0.000 


0.000 


-16.406 


252.594 


Table  18.  NUMERICAL  DATA  CORRESPONDING  TO  CASE  ALP5-C 


Case  5c:  Three  targets  located  at  random  positions  -  targets  one  and  two  have  a 
common  spectral  line  at  6000  Hz  -  six  harmonics  present  in  output  electrical  signal  - 
ALP  algorithm  is  used  -  1000  iterations  -  no  noise 


Average  Estimation  Error 
(des) 


Location  (deg) 

00 

77 

307 

5 

137 

49 

38 

77 

307 

5 

137 

5 

137 

49 

38 

0.094 


-0.038 


0.094 


0.094 


0.006 


2.615 


0.498 


0.012 


-0.008 


0.012 


0.012 


-0.328 


-3.84 


2.516 


i 
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